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ABSTRACT 

The linear polarisation of the Crab pulsar and its close environment was derived from 
observations with the high-speed photo-polarimeter OPTIMA at the 2.56-m Nordic 
Optical Telescope in the optical spectral range (400-750 nm). Time resolution as short 
as 11 /is, which corresponds to a phase interval of 1/3000 of the pulsar rotation, and 
high statistics allow the derivation of polarisation details never achieved before. The 
degree of optical polarisation and the position angle correlate in surprising details with 
the light curves at optical wavelengths and at radio frequencies of 610 and 1400 MHz. 
Our observations show that there exists a subtle connection between presumed non- 
coherent (optical) and coherent (radio) emissions. This finding supports previously 
detected correlations between the optical intensity of the Crab and the occurrence of 
giant radio pulses. Interpretation of our observations require more elaborate theoretical 
models than those currently available in the literature. 

Key words: radiation mechanisms: non-thermal - pulsars: general - pulsar- 
individual: the Crab pulsar - techniques: polarimetric - instrumentation: polarimeters 



1 INTRODUCTION 

The Crab pulsar and its pulsar wind nebula (PWN) are two 
of the most intensively studied objects in the sky. The com- 
pact remnant of SN1054, a cornerstone of high energy astro- 
physics, is one of the youngest and most energetic pulsar and 
its pulsed emission has now been detected th roughout the 
electromagnetic spectrum from about 10 MHz (|Bridlelll970h 
up to > 25 GeV wit h some evidence (3.4a) of pulsed emis- 
sion above 60 GeV (|Aliu et all 120081 ). The PWN was mea- 
sured even up to e nergies of ~ 100 Te V l|Aharonian et all 
|2004 lAllenl 120071 : lAllen et all 120071 ). The pulsar and its 
nebula are predominantly sources of non-thermal radiation 
(synchrotron, curvature and inverse Compton processes), 
which is indicated not only by the broad-band spectral con- 
tinua but also by the strong polarisation of these emissions. 
The outstanding brightness of the Crab across the electro- 
magnetic spectrum makes it the ideal target to investigate 
polarisation at all wavelengths where suitable instrumenta- 
tion is available. 

Before we focus on the Crab we want to introduce 
the subject of polarisation in pulsars by a short review 



of the total population: most pulsars were discovered in 
the radio regime (presently about two thousand objects 
are listed in the Australia Telescope N ational Facility Pul- 
sar CatalogueFl. iManchester et al.ll200o] ) and for for a large 
fraction of them po l arisation studies were performe d (e.g. 
iGould fc Lv"iie1ll998l : lK"arastergiou fc Johnston! 120061 ). Most 
radio pulsars were found to show strong linear polarisation, 
including often a characteristic a swing of the position angle 
(P.A.) in an S-like shape near the pulse centre. This swing is 
interpreted in the 'rotating ve ctor model' (hereafter: RVM, 
iRadhakrishnan fc Cooke] 1 19691 ) as a projection of the mag- 
netic field line at the point of emission onto a plane perpen- 
dicular to the observer's sight-line. The point of emission is 
usually assumed to be in the polar cap region of the pulsar 
where a regular dipolar field line points with a small an- 
gle (beam width) towards the observer. The free parameters 
of this simple geometrical model are the inclination angle 
between the axes of rotation and magnetic dipole, and the 
viewing angle between the line of sight and the rotation axis. 
Analyses of radio polarisation from ma ny pulsars (e.g. ref- 
erences m iLvne fc Graham-Smith] 120061 ) showed that most 
sources can be described in a polar cap model with an RVM 
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Table 1. High energy polarisation measurements of pulsars and their nebulae. Polarisation degree is given as a percentage, and position 
angle in degrees (North, 0° to East, 90°). 



Band 




Pulsar / Polarisation 


Nebula (near PSR) 


Ref. 


Optical (V~ 


' 16.6) 


Crab 

phase-resolved (Fig. 
phase-averaged 9.8% ± 0.1% and 109.5° ± 0.2° 


9.7% ±0.1% ; 139.8° ±0.2° 
(< 5" from PSR) 


[1], [2], [3] 


UV 




phased-resolved, similar to optical - 
r.L). and r\A. [Mr ana lr) 




W 


X-ray 




only upper limits 


19.2% ± 0.9%; 155.8° ± 1.4° 


[5], [6] 


Hard X-ray/ soft 7-ray 


off-pulse (phase: 0.52-0.88) > 72%; 120.6° ± 8.5° 
phase-averaged 47%+^°; 100° ± 11° 




[7] 
[7] 


7-ray 




off-pulse (phase: 0.5-0.8) 46% ± 10%; 123° ± 11° 




[8] 


Optical (V~ 


' 22.5) 


B0540-69 

phase-averaged: ~ 5%, no error quoted; phase-resolved: < 15% 


5.6% ± 1.0% ; 79° ± 5° 


[9], [10], [11] 


Optical (V~ 


' 23.6) 


Vela 

phase averaged: 9.4% ± 4%, 146° ± 11° 




[9], [12] 


Optical (V~ 


' 25) 


B0656+14 

double peak light curve; P.D. bridge ~ 100%, peaks ~ 0% 
P.A. sweeps in agreement with RVM 




[13] 


Optical (V~ 


- 25.7) 


B1509-58 

phase-averaged: ~ 10.4% (very uncertain, no error quoted) 




[9] 



fll ISmith et alj lll98Sft . [21 ISlowikowska et alj ll2008t). [31 this pap er, [ 41 iGraham-Smith et all |l996l l . [51 ISilver et al.l (Il978h. [ 61 
IWeisskopf et al.l lll978l) , [71 iForot et al.l 1I2OO8II. [ 81 bean et all d2008h. [9l lWagner fc SeifertJ <200ch . [101 iMiddleditch et al.l lll987fl , [11] 
IChanan fc Helfandl <199C|) . [121 iMienani et alj 1120071) . [13] iKern etafl ( 12003^ 



and the intrinsic polarisation could be as high as 100%. Sev- 
eral pulsars however, especially at radio frequencies above 
several GHz, show reduced or complex polarisation signa- 
tures with indications of depolarisation or switching between 
modes of polarisation. These complications could arise from 
the subtle and often highly variable and unstable nature of 
the coherent processes underlying the radio emission. 

The situation changes when going to higher energies, 
i.e. optical, X- and 7-ray regimes. Incoherent single parti- 
cle radiation processes in the magnetosphere provide the 
high-energy pulsar emission (thermal emission from the neu- 
tron star surface is not considered here). Precise polarisation 
measurements of polarisation degree (P.D. jf| and P.A. as a 
function of pulsar phase in a wide range of energies should 
provide deep insight into the pulsars emission mechanisms, 
particle spectra, and emission site topologies. Fast X- and 7- 
ray polarimetry from space borne instruments is presently of 
very limited sensitivity. Therefore results have been reported 
only for the brightest pulsar and its PWN, i.e. the Crab. In 
the optical domain the situation is somewhat reversed: here 
the polarimeters are quite sensitive but the optical magni- 
tudes of most pulsars are very faint. Therefore, despite the 
increasing number of optical pulsars (fourteen are known), 
only for five of them attempts to measure the pulsar and/or 
nebular polarisation were made. Table Q] summarises these 



2 The terms: polarisation degree, P.D. and p as well as position 
angle, P.A. and 8 are used interchangeably. 



efforts. For three pulsars phase-resolved measurements were 
performed; for the remaining two only phase-averaged re- 
sults are available. Again the Crab, being the brightest ob- 
ject, is the exception with a multitude of optical polarisation 
measurements. 

The scarcity of optical polarisation measurements, ei- 
ther in a phase averaged mode or with time resolution 
of the pulsar rotation, is clearly due to the large differ- 
ence in magnitude going from the brightest pulsar, i.e. 
the Crab, to the next, B0540-69, which is ~ 6 magni- 
tudes fainter. Therefore only for the Crab pulsar fully 
phase resolved polarisation measurements have been pos- 
sible so far. The first optical phase-resolv ed linear polar- 
isation observations of the Crab pulsar ([Wampler et all 
1 19691 ; ICocke et all Il97fj ; iKristian et all Il970l ) showed that 
the polarisation angle sweeps through each peak and the 
polarisation degree decreases and then increases within 
each pulse, reaching the minimum shortly after the pulse 
peak. These early observations were limited to the main 
pulse (MP) and inter pulse (IP) phase ranges only. Inter- 
pretation of the polarisation pattern behaviour has been 
made in terms of geometrical models ([Wampler et al.lll969l; 
[ Radhakrishnan fc Cooke! 19691 ; ICocke et al.l 19731 ; Ferguson! 
Il973l ; iFerguson et al.lll974T ). some of which offer the possi- 
bility of determinin g the magnetospheric location of the ra - 
diation source (e.g. ICocke et al.lll973l ; IFerguson et~aLlll974l ) . 

For a long time it was thought that optical radiation 
of the Crab pulsar persisted only through both peaks, the 
MP and IP, and in the bridge region between them, but not 



in the phase range following the IP and pre ceding the MP. 
Several phase resolved imaging observations ([Peterson et al.l 
1 19781 : IPercival et aDll993l : Icolden et al.ll2000h showed how- 
ever that radiation persists throughout the whole pulsar ro- 
tation. They found that the minimum intensity, consistent 
with coming from an unresolved source at the pulsar posi- 
tion, occurs immediately before the start of the MP within 
a phase range of approximately 0.779 and 0.845. This level, 
often called the DQ_| level of the Crab pulsar, has been vari- 
ously determi ned with an intensity relative to the MP maxi- 
mum of 3.6% (iPeterson et al.lll97Sft < 0.9% (|Percival et all 
1 1993! ). or ~ 1% (|Golden et al.ll2000D . 

After this discovery, the linear polarisation of optical ra- 
diatio n fro m the Crab pulsar was measured bv lJones et al.l 
<|l98lf ) and ISmith et all (|l988h . Both results confirmed the 
previous observations and additionally showed the polar- 
isation characteristics during both off-pulse phase ranges 
(bridge and DC). However, for the DC phase range the re- 
sults were not conclusive. The polarisation degree was at a 
level of 70% an d 47 ± 10% and the positi on angle was ~ 118 ° 
and ~ 130° for I Jones et ail l|l98ll ) and ISmith etafl l|l98Sh . 
resp ectively. 

IWeisskopf et al.l l| 19781 ) measured the linear polarisation 
of the X-ray flux from the Crab nebula at energies of 2.6 keV 
and 5.2 keV, independent of any contribution from the pul- 
sar. Within a field of view of 3° the X-rays are polarised at 
a level of 19.2% - 19.5% and a P.A. of 156° - 1 52°. For the 
pulsa r itself only upper limits could be derived l|Silver et al.l 
1 19781 ). No evidence for an energy dependence of the X- 
ray polarisation was found. Close to the pulsar the nebu- 
lar optical polarisation is quite uniform (at ~ 9% — 11%) 
but the position angles change steadily with radial distance; 
2" — 3" from the pulsar the mean value is around 140° but 
beyond 5" the posi tion angle exceeds 155° and it is very 
location-dependent (|McLean et alj 1983). The X-ray results 
for the nebula are in good agreement with the optical po- 
larisation measurements which yield a 19% polarisation at 
162° for the central regio n of the nebula within 0'.486 radius 
l|Oort fc Walravenlll956h . The similarity of the optical and 
X-ray results indicates that the polarisation is independent 
of energy over a broad spectral range. Recently 7-ray phase- 
averaged (off-pulse phase region, i.e. 0.5-0.8) polarisation 
m easurements of the Crab nebula and pulsar were reported 
bv lDean et ail (|2008l ) for the energy range 100 kV- 1 MeV. 
The P.D. was found at 46%±10% and the P.A. is 123°±11°, 
which indicates that this energetic radiation is dominated by 
emission from the inner nebula arou nd the pulsa r . This re- 
sults has been recently confirmed bv lForot et al.1 (|2008h by 
using the IBIS instrument on board of the INTEGRAL tele- 
scope to obtain the linearly polarised emission for energies 
between 200 — 800 keV and extend polarisation characteris- 
tics for pulses and bridge phases (Tab.[T]). 

We structure the paper as follows: firstly we describe our 
instrumentation (§2). In §3 and §4 the observations, data re- 
duction methodology and measurements of polarimetric and 
photometric standards are described. In §5 we present the 
results of the Crab nebula and its pulsar. Also in §5 we 



3 Direct Current or continuous current, in this case refers to a 
continuous emission component that could be present throughout 
the whole rotational phase. 
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Fiber Bundle on HST image (scaled to NOT) 
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Figure 1. An enlargement of the HST image of the inner Crab 
nebula (source R. Romani, priv. comm.). North is up; East is left. 
The pulsar is identified with the lower / right of the two stars near 
the geometric centre of the nebula. A small arc-like feature (the 
inner knot) which is clearly resolved in the HST image is located 
(y/65 to the SE of the pulsar. The OPTIMA fibre bundle, centred 
on the pulsar and scaled with the NOT focal plane scale, is over 
plotted. 



show the pulsar polarisation characteristics after subtrac- 
tion of an unpulsed component and discuss time alignment 
between optical and radio wavelengths. In §6 we present a 
discussion of theoretical models and compare them to our 
measurements and end with a summary and conclusions. 



2 INSTRUMENTATION 

Our goal was to obtain polarisation characteristics of the 
Crab pulsar with very high time resolution as a function of 
the pulsar rotational phase. We used the 2.56-m Nordic Op- 
tical Telescop43 (NOT) for these measurements. This facility 
not only provides a mirror with reasonably large collecting 
area but also has an excellent performance in offset guid- 
ance and telescope control, which is important for instru- 
ments with very small entrance apertures. Moreover, con- 
trary to many larger optical telescopes, guest instruments 
can be used at NOT. For the Crab pulsar observations car- 
ried out during November, 23-27 2003, we used the high- 
speed photo-polarimeter OPTIMA0 (Optical Pulsar TIMing 
Analyser). We shortly describe the system below. An ex- 
te nded description of the OPT IMA instrument is presented 
bv lKanbach et al.l (|2003l . 120081 ). 



2.1 OPTIMA 

The basic requirement for a high-speed photometer is to se- 
lect and time tag individual photons from a celestial source 
with high efficiency and precision. In the case of OPTIMA, 
which was primarily designed for studying the optical light 
curves of faint pulsars and other highly variable targets, the 



4 http : //www. not . iac . es/telescope/technical-details .html 
http : //www.mpe .mpg.de/gamma/ instruments /optima/ 



4 A. Slowikowska et at 



source flux is isolated in the focal plane of the telescope by 
the use of an optical fibre pick-up which acts as a diaphragm. 
The target fibre is central to a hexagonal bundle of identi- 
cal fibres which measure the sky background and, in the 
case of the Crab nebula, the close environment of the pulsar 
(Fig. [T}. OPTIMA uses commercial single photon count- 
ing modules with Avalanche Photodiodes (APDsjfj. These 
APD counters have quantum efficiencies peaking at 70% for 
A ~ 700 nm and a wide range with Q.E. > 20% from 440 to 
980 nm. 

The light from the telescope at the Cassegrain focus 
is incident on a slant mirror with an embedded bundle of 
optical fibres. Optionally one can insert filters or a rotating 
linear polariser into the incoming beam. The field around the 
fibres, visible in the mirror (typical size 2' x 3'), is imaged 
with a target acquisition camera (type AP6, from Apogee 
Instruments). We use the following nomenclature for the 
fibre channels: channel for the central fibre, channels 1 
to 6 for the ring fibres and channel 7 for a sky background 
fibre, located about 1' off-set from the target fibre. For the 
observations we used tapered fibres with 320 /im diameter at 
the pick-up and an exit diameter at the detector of 100 fim. 
The single fibre size in the focal plane is equivalent to a 2'.' 344 
resolution at the 2.56-m NOT telescope. Although the pure 
silica core fibres transmit light from 400 nm to ~950 nm 
with better than 99%, the net transmission of the tapered 
fibres is lower and appears to range between 60% and 90%. 

The timing of individual photons is controlled by signals 
from the Global Positioning System (GPS) to an absolute 
accuracy of ~ 2 [is, although the readout system limits the 
resolution to ~ 4 /is. The OPTIMA detector is operated 
with two PCs and is autonomous except for the need to 
have a good telescope guiding system. 

2.2 Rotating Polarisation Filter 

For the present observations OPTIMA was equipped with 
a rotating polarisation filter (RPF) in the incoming beam 
above the fibre pick-up so that all fibre channels and the 
CCD image are fully covered. The polarising filter (Type 
10K by Spindler & Hoyer) is mounted on a precision roller 
bearing and is rotated with typical frequencies of a few Hz 
(the averaged frequency over the whole observations was 
3.4 Hz). Incoming linearly polarised light is then modulated 
at twice the rotation frequency of the filter. The reference 
position of the filter is given by a signal from a magnetic 
switch (Hall sensor) which is registered and timed in the 
same way as a photon event and stored in a separate DAQ 
channel. The position of the polarising filter for any photon 
event is then derived by interpolating the photon arrival 
time between the preceding and the following Hall sensor 
signal, i.e.: 

a RPF (TOA) = tTOA - h x 360°, 

ta — tb 

where okpf(TOA) is the RPF angle at which the event was 
observed, troA is the event time of arrival, and tb and t a 
are the recorded times of the Hall signal sensor before and 
after troA- Since, during one turn of the RPF all possible 

6 Perkin-Elmer SPCM-AQR-15FC 



Table 2. Fibre flat field correction factors 



Channel 


Factor 





1.000 ±0.012 


1 


1.009 ±0.013 


2 


1.074 ±0.013 


3 


0.899 ± 0.012 


1 


0.850 ±0.011 


5 


0.956 ±0.012 


6 


1.140 ±0.014 



polarisation angles are measured two times, for all values 
bigger than 180° exactly 180° was subtracted. 

Slight irregularities in the rotation frequency of the 
RPF that occur on time-scales longer than fractions of a 
second, e.g. due to supply voltage drifts or mechanical re- 
sistance changes in the bearing and motor, can thus be cor- 
rected with sufficient accuracy. The RPF was tested in the 
lab with unpolarised and linearly polarised light to ensure 
and prove that the OPTIMA fibres and detect ors have no in - 
trinsic systematic response to polarised light (|Kellnerll200"2l ). 

The polarising filter modulates the incoming light ef- 
fectively only over a wavelength range of about 470-750 nm. 
Since the APD response (QE> 20%) extends from about 
450 nm to 970 nm and no wavelength information of the in- 
dividual recorded events is available, it is necessary to block 
radiation outside the filter modulation range. Such photons, 
especially towards the near IR, are not modulated and would 
decrease the estimate for the degree of polarisation. There- 
fore, we inserted an IR blocking filter that cuts the wave- 
length range at about 750 nm. 

The method to derive the pol arisation charac t eristic s 
follows the approach described by ISparks fc Axon! (|l999h . 
Details are given in Appendix|X]to this paper. 

3 OBSERVATIONS 

The Crab observations were performed on November 23-27 
2003 (52966 - 52971 MJD). They consists of 160 pointings 
of ten minutes each. About one third of the total number 
of counts falls in the central fibre, whereas around 5 to 10 
percent are registered in each of the background fibre (chan- 
nels marked from 1 to 6 in the Fig. [TJ. The count rate for 
the central fibre decreases when the seeing deteriorates and 
photons spill out into the ring channels (1 ± 6). In order 
to screen for observation intervals of acceptable seeing we 
therefore apply a cut on the fraction of total counts in chan- 
nel at the level of 30%. After screening for good seeing and 
proper pointing, our data base resulted in 83 files, equivalent 
to 13 hours and 43 minutes of total exposure. 

4 DATA REDUCTION 
4.1 Flat field correction 

For flat field correction factors of all channels we used 
two sets of dark sky observations, i.e. pointing at a = 
00 h 07 m 02 s 00j s = + 73°03'28'.'00. Both observations were 
taken on November 25. The data acquisition started at 
21:39:10 UTC and 21:49:13 UTC. The exposure time 
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Table 3. HST polarisation and pho tometric standards chosen 
from the list of iTurnshek et ail l|l990h 



Name / Coordinates FK5 

Spectral Type 



Comments 



HST Polarisation Standards 

BD+64 106 00 h 57 m 36?70 V = 10.34 

B1V +64° 51' 34'.'9 p = 5.65% ± 0.053% 

9 = 96?8 

G191B2B 05 h 05 m 30?61 V = 11.79 

WD +52° 49' 51'.'9 p = 0.09% ± 0.048% 

HST Photometric Standard 



GD50 

WD 



03 h 48 m 50? 20 
-00° 58' 31'.'2 



V = 14.06 



amounted to 600 and 300 seconds for the first and second 
observation, respectively. The detected count rates behind 
the filters fluctuated from 180 Hz for channel 4 up to 240 Hz 
for channel 6. To obtain the flat field correction factors we 
assumed that this value was 1.0 for the central fibre, and we 
scaled the count rates of the other fibres to the central one. 
We binned the data in ten seconds intervals, and then cal- 
culated the average and standard deviation for each single 
channel of the ring fibres. The resulting numbers are shown 
in Tab. [2] 



4.2 The HST polarisation standards 

We performed observations of stars with well known polari- 
sation in order to understand the intrinsic polarisation and 
the response of the instrument. Since OPTIMA is very sen- 
sitive and limited in its capacity for data acquisition at high 
rates (count rates above 37 kHz lead to n oticeable, but cor- 
rectable, pile-up effects, [Miihlcgger 20 0if) we selected thre e 
of the weakest stars from the list of Turnshek et al.l {l990); 
two polarisation and one photometric standards (Tab. [3}. 
The first star of the two polarisation standards is highly 
polarised, on the level of 6%, whereas the second one has a 
very low polarisation degree of ~1 %o- Unfortunately for our 
measurements both of them are quite bright by OPTIMA 
standards. Therefore, we additionally performed optical po- 
larisation mea surements of a d i mmer star, i.e. a photometric 
standard from lTurnshek et al.l (|l990T ) list, expecting its light 
not to be polarised. 

We observed the HST polarisation standard BD+64 106 
for 270 and 350 seconds on the 25th and 27th of November, 
respectively (Tab. [4}. The observing conditions during the 
two exposures were very different, requiring a detailed dis- 
cussion of the results. On Nov 25, 22:06 - 22:12 UTC the 
average seeing was on the level of l'.'27 (RoboDIMM mea- 
surements), whereas for the time span 21:33 - 21:40 UTC 
on Nov 27 no information about the seeing conditions was 
available from the RoboDIMM telescope. At the beginning 
of this night the weather was very good. But later, during 
less than two hours, between 20:00 and 22:00 hours, the hu- 
midity changed from 20% to almost 60%, and the wind speed 
increased up to 15 m/s. Therefore, the seeing on Nov 27 was 



probably worse than during Nov 25, likely more than 2'.'0. 
This difference in seeing causes different count rates in the 
central and ring fibres. On 25 Nov the count rate in channel 
was about 190 kHz, and in the ring fibres ss 20 kHz, while 
during Nov 27 we had 150 kHz (centre) and w 50 kHz (ring). 
As a consequence of the high count rates (3> 37 kHz) the 
central channel was severely affected by pile-up in these ob- 
servations. For the second set of data (27 Nov) it happened 
that some of the ring fibres, especially channels 1 and 6, were 
also saturated. The values of the polarisation degree and the 
position angle for all channels and for all sets of observations 
are gathered in Tab. [3] It is important to notice that, even if 
the central fibre was suffering pile-up, the calculated position 
angle for each single observation of BD+64 104 is very much 
constant. This is caused by the fact that the position angle, 
being the phase of the maximum of the modulated incoming 
light, is not strongly dependent on the pile-up effect. Sat- 
uration is however important when the amplitude of RPF 
modulation, corresponding to the degree of polarisation, is 
to be measured. We can observe this by comparing results 
in Tab. 2] For the central fibre there is significant difference 
in the polarisation degree for the first and second data set. 
When the dispersion of the stellar image due to seeing was 
larger the star light was more smoothly distributed among 
all channels, therefore an increase of p in the central fibre is 
observed. Seeing conditions, being equivalent to the count 
rates in the ring, did not affect significantly the measured 
values of p in the ring fibres. Only in the case of channel 1 
the polarisation degree is much smaller during the Nov 27 
observations than Nov 25. This is strongly connected with 
the fact that, during the second observation, this channel, 
among all of the ring channels, was the most saturated one. 

Taking into account the resulting 6 and p (Tab. 2} of 
the HST polarisation standard BD+64 104 (being a very 
bright star as for the OPTIMA instrument) we can conclude 
that for the purpose of calibrating the north direction of our 
instrument it is safe and valid to use values of the polarisa- 
tion angle from the central fibre (Tab. [4] bold faced text). 
The averaged value of 6 i bd+64104 amounts to 96.74° ±0.37°. 
This result was obtained by shifting the intrinsic RPF an- 
gles (the position of the Hall sensor indicates the origin of 
the intrinsic angles) by 92°. The angles 9 given here and 
hereafter in the text and tables are the E— vector position 
angles relative to celestial north (N to E). Due to the pile- 
up in the central fibre the best estimate of the polarisation 
degree of BD+64 104 comes from the sum of the light de- 
tected in the ring fibres. The obtained value, 4.72% ±0.30%, 
is som ewhat smaller than the value given by I Turnshek et ail 
l|l990h . 5.65%±0.053%, but one should take into account the 
systematic problems encountered with such a bright stan- 
dard star. In this case the bad seeing (high value) is an ad- 
vantage, because it naturally defocuses our target. We also 
performed observations of the unpolarised HST standard 
G191B2B (Tab. 13 and [SJ. Data for this target were taken on 
26 Nov, 21:33 UTC for about 10 minutes with seeing being 
on the level of l'.'O. We found a very small polarisation de- 
gree for this star, as expected, and thus the position angle 
has little meaning. During these observations the count rate 
was on the level of 150 kHz and <40 kHz in the central 
and ring fibres, respectively. We should point out that chan- 
nel was affected by pile-up and the result 0.05% ± 0.02% 
(Tab. [5]l might be biased, even if it is in good agreement with 
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Table 4. Measurements of the HST polarisation standard BD+64 104. Date and time of observations, as well as exposure length in 
seconds are given in the first, second and third row, respectively. Obtained values of P.D., p, (upper value) and P.A. (9, bottom value) are 
given for each OPTIMA detector channel from to 6 (p and 9 are in percent and degrees, respectively). Additionally, last row (assigned 
as '1 4- 6' ) gives p and 9 values obtained after averaging the Stokes parameters from all ring channels. Mean gives p and 9, as well as 
their standard deviations, after averaging all existing data of BD+64 104. The obtained position angle of BD+64 104 (bold faced text 
in the row of channel 0) was later used for the OPTIMA instrument calibration. 



BD+64 104 


Obs. Date 


25 Nov 2003 


27 Nov 2003 




Mean 


lime [U 1 C] 


99-07-01 


99-OQ-zlfi 


91 -^-^ 


91 -^fvOS 


21:38:34 




Expo, [s] 


140 


130 


110 


120 


120 




Channel 

















1.30 


1.35 


1.98 


1.86 


1.67 


1.63 + 0.30 




96.6 


96.8 


96.2 


96.9 


97.2 


96.74+0.37 


1 


4.57 


4.64 


3.71 


3.60 


4.21 


4.15 + 0.48 




100.3 


96.6 


98.0 


96.5 


98.3 


98.14 + 1.62 


2 


5.16 


5.17 


4.59 


4.84 


5.57 


5.07 + 0.37 




96.0 


94.7 


98.0 


96.8 


96.8 


96.46 + 1.22 


3 


5.38 


5.16 


5.49 


5.08 


5.04 


5.23 + 0.20 




94.3 


92.9 


95.5 


95.9 


97.4 


95.20 + 1.70 


4 


4.67 


4.34 


4.14 


4.98 


4.27 


4.48 ± 0.34 




99.1 


98.6 


98.5 


100.5 


103.5 


100.04 ± 2.09 


5 


5.51 


5.26 


5.53 


5.58 


5.38 


5.45 + 0.13 




96.7 


94.3 


97.9 


96.2 


97.1 


96.44 + 1.35 


6 


4.62 


5.00 


5.06 


4.52 


4.80 


4.80 + 0.23 




98.6 


98.7 


93.6 


93.6 


91.7 


95.24 + 3.21 


1H-6 


5.06 


4.95 


4.47 


4.37 


4.73 


4.72 + 0.30 




96.7 


95.1 


96.8 


95.8 


96.4 


96.16 + 0.71 



Table 5. Measurements of the HST polarisation standard G191B2B and the HST photometric standard GD50. Date and time of 
observations, as well as exposure length in seconds are given in the first, second and third row, respectively. Obtained values of P.D., p, 
arc given in percent for each OPTIMA detector channel from to 6. Additionally, last row (assigned as '1 4- 6') gives p values obtained 
after averaging the Stokes parameters from all ring channels. Mean gives p and standard deviations after averaging all existing data of 
G191B2B from Nov 26th. In the last column the results of GD50 measurements are given. 









G191B2B 






GD50 


Obs. Date 




26 Nov 2003 




Mean 


27 Nov 2003 


Time [UTC] 


21:33:08 


21:36:22 21:41:01 


21:43:53 




21:52:17 


Expo, [s] 


160 


150 


130 


130 




590 


Channel 

















0.04 


0.07 


0.03 


0.06 


0.05 ± 0.02 


0.08 


1 


0.37 


0.48 


0.47 


0.21 


0.38 ± 0.13 


0.30 


2 


0.75 


0.41 


0.65 


0.85 


0.67 ± 0.19 


0.44 


3 


0.84 


0.46 


0.80 


0.72 


0.71 ± 0.17 


0.61 


4 


0.30 


0.66 


0.69 


0.47 


0.53 ± 0.18 


0.62 


5 


0.16 


0.56 


0.36 


0.16 


0.31 ± 0.19 


0.23 


6 


0.14 


0.21 


0.04 


0.42 


0.20 ± 0.16 


0.53 


1^6 


0.35 


0.34 


0.30 


0.32 


0.33 ± 0.02 


0.22 
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Table 6. The Crab pulsar radio ephemeris 



Parameter 



Value 



R.A. (J2000), «2000 
DEC (J2000), 5 2 ooo 
Valid range (MJD) 
Epoch, t (TDB MJD) 
vo (Hz) 

u (1CT 10 Hz s" 1 ) 
u (10- 20 Hz s" 2 ) 



05 ?l 34 m 31?972 
+22°00'52'/07 
52944-52975 
52960.000000296 
29.8003951530036 
-3.73414 
1.18 



the value 0.09% ± 0.048% given by iTurnshek et all (|l990h . 
The averaged value of P.D. from the background fibres is 
higher. It amounts to 0.33% ±0.02%, and is probably partly 
contaminated by the polarised sky background. 

For the pu rpose of the ins t rume nt calibration we also 
chose from the ITurnshek et al.1 (Il990l ) atlas of HST photo- 
metric, spectroscopic, and polarimetric calibration objects 
a star with brightness that should assure count rates lower 
that the pile-up threshold, i.e. GD50 (Tab. [3}. It is a photo- 
metric standard, thus we expect its light not to be polarised. 
Observations of this target were performed on Nov 27th at 
21:52 UTC with recorded count rates in the central chan- 
nel of 10-15 kHz and in the ring channels < 1 kHz. Results 
are shown in the Tab. [5] The source had polarisation degree 
below l%o (measured in the central fibre), whereas the mea- 
sured p within the ring fibres was about three time greater. 
These values correspond to the sky background polarisation. 
However, the sky background could have been more highly 
polarised. The contribution of the night sky to the observed 
count rates was on the level of one quarter in each fibre 
only, i.e. the background ring fibres were contaminated by 
the source. Thus, the signal measured in the ring fibres was 
most probably depolarised and could have higher polarisa- 
tion degree in reality. 



4.3 Raw data binning for the Crab analysis 

For each incoming and detected photon the OPTIMA data 
acquisition system stores its Time Of Arrival (TOA) in a 
compact proprietary binary format with a resolution of 4 /is . 
By using the OPTIMA system software (|Straubmeierll200ll ) 
one can obtain TOAs in units of Julian Date (JD). In or- 
der to apply the rotational model of the pulsar derived from 
radio observations to our optical photons, it is required to 
measure the TOAs in an 'inertial observer frame'. Using 
the NOT position coordinates, the Crab pulsar coordinates 
(Tab.[5| and as plane tary ephemeris the JPL DE200 tabula- 
tions i|standishl fl982) we transformed the recorded TOAs to 
the commonly used, and (for these purposes) inertial frame, 
i.e. to TOAs at the solar system barycentre. 

Parameters of the rotational model of the Crab pulsar 
(Tab. [6} obtained from the radio observations are regularly, 
once per month, published by the Jodrell Bank Observa- 
tory pulsar group in the Crab Pulsar Monthly Ephemerifl 



tory pulsar group 
|Lvne et al.lll993l ) 



For each recorded event (TOA) the cor- 
responding pulsar phase and the phase of the RPF is cal- 
culated. We sort the individual events into a 3D data ar- 



http: //www. jb .man. ac .uk/~pulsar/crab. html 
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Figure 2. The Crab pulsar light curve obtained from photons 
recorded by the OPTIMA central fibre APD (channel 0). The 
components of the light curve are indicated as follows: main pulse 
MP, inter pulse IP, non-zero intensity level between MP and IP, 
i.e. bridge, as well as the DC region, previously known as the 
'off-pulse' component in in the phase range 0.7729 — 0.8446 . The 
upper panel shows the raw light curve including counts from the 
nebula and the sky background. In the lower panel the lightcurve 
is shown in a logarithmic scale after subtraction of the nebular 
and sky backgrounds. This display shows clearly the DC intensity 
level. Two rotation periods are shown for clarity. 



ray. The first dimension is given by the rotational phase of 
the pulsar (binned with various resolutions), the second one 
by the phase of the polarisation filter at each photon TOA 
(binned in 1° intervals), and the third dimension is the fibre 
number or number of DAQ channels (seven channels). 

To obtain an 'unpolarised' light curve, e.g. with a time 
resolution of 67 /is which assures a good S /N ratio even for 
the phase ranges where the intensity is very low, the data 
array is constructed with 500 bins per rotational period of 
the Crab, and the events spread out over the RPF phases 
are all added up. The Crab light curve, in its raw form as 
the measured count rate in the central fibre without back- 
ground subtraction, is shown in the top panel of Fig. [2] We 
indicate the components main pulse (MP), inter pulse (IP), 
non-zero intensity level between two peaks - bridge, and the 
DC region, previously known as the so called 'off-pulse' com- 
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Table 7. The Crab nebula's P.D. and P.A. at the DC phase range 





Channel 


p[%] 


e n 





33.08(25) 


118.8(0.2) 


1 


9.05(8) 


141.3(2) 


2 


10.30(8) 


142.5(2) 


3 


8.77(7) 


148.4(2) 


1 


9.35(8) 


136.6(2) 


5 


11.50(8) 


133.1(2) 


6 


10.26(7) 


137.9(2) 


l-r-6 


9.71(8) 


139.8(2) 



ponent. For our further analysis we define the DC compo- 
nent as the counts between phases 0.7729 — 0.846 , whic h 
is in accordance with the findings of IPercival et al.l l|l993h . 
The same light curve after background subtraction (method 
described in the next section) is shown in the lower panel of 
Fig. [2] with a logarithmic scaling to enhance the visibility of 
the DC level. The higher S/N ratio in the peaks allows us 
to use better time resolution for these phase ranges, i.e. up 
to 3000 bins per period, which corresponds to a bin of 11 /is 
interval (Figs. U [5] |6j). 

5 RESULTS FOR THE CRAB NEBULA AND 
PULSAR 

5.1 Nebular contribution 

The Crab synchrotron nebula is a relativistic magnetised 
plasma that is powered by the spin-down energy of the 
pulsar. It was the first recognised astronomical source of 
synchrotron radiation. The synchrotron nature of the ra- 
diation was co nfirmed by optical polarisation observations 
l|Woltienll957h . The conversion efficiency of the nebula is 
quite high, with 10%-20% of the spin-down energy released 
by the pulsar appearing as synchrotron radiation. The in- 
ner synchrotron nebula is a region consisting of jets, a torus 
of X-ray emission, small-scale variations in polarisation and 
spectral index, and complexes of sharp wisps. Most theo- 
retical models associate the sharp wisps seen at visible and 
radio wavelengths with the location of the shock wave be- 
tween the pulsar and the synchrotron nebula. Closer to the 
outer boundary of the nebula are the filaments - the chemi- 
cally enriched material ejected during the supernova explo- 
sion observed by Chinese astronomers in 1054. 

For a long time the imaging investigations of the Crab 
nebula were limited in a fundamental way by the spatial 
resolution of the detectors, not being able to reveal a wealth 
of subarcseconds structures. A breakthrough in the optical 
studies of the structure of the Crab nebula was undertaken 
by using the Wide Field and Planetary Camera 2 (WFPC2) 
on board of the Hubble Space Telescope. iHester et al.l (|l995l ) 
performed this observation with O'.'l resolution. Below we 
shortly describe the most important HST results that are 
important in the context of this paper. 

They discovered a bright knot of visible emission lo- 
cated 0'.'65 to the south-east of the pulsar, along the axis 
of the system (Fig. [TJ. This inner knot, along with a sec- 
ond similarly sharp but fainter knot (hereafter outer knot) 
located at a distance of 3" 8 from the pulsar, lies at an ap- 




Figure 3. Polarisation of the Crab nebula as measured by each 
single fibre at the minimum phase of the pulsar light curve. It is 
graphic representation of results given in Tab. [7] The axes are cali- 
brated in arc seconds centred on the pulsar. The scal ed OPTIMA 
fibre bundle is over plotted on the figure given in ISmith et al.l 
Il988l . The aperture of a single fibre is 2'.'35. The pink line for cen- 
tral fibre and the red lines for the ring fibres represent our mea- 
surements with the rotating pola roid filter. The black lines are 
the results that ISmith et al,lll988l determined with a polarimeter 
based on a rotating half-wave plate and a Foster prism. 



proximate position angle of ~ 115° east to north. Both knots 
are aligned with the X-ray and optical jet to the south-east 
of the pulsar and are elongated in the dimension roughly 
perpendicular to the jet direction, with lengths of about a 
half arc second. Both, the inner and outer knot appear to 
be present but not well resolved in the images of the Crab 
nebula previously taken by ground based telescopes. Fig. [T] 
is an enlargement of the co-added HST WFPC2 images of 
12 observations of the Crab nebula in between 2000 and 
2001 with the F547 M filter (observations group index 6 of 
iNg fc Roman! 2006, it was kindly supplied by Roger Ro- 
mani, private communication). The pulsar is identified with 
the lower/right of the two stars near the geometric centre of 
the nebula. The OPTIMA fibre bundle centred on the pulsar 
and scaled with the NOT focal plane scale is over plotted. 
It is clear that we are not able to resolve the Crab pul- 
sar from the inner knot within the central fibre. Moreover, 
with seeing being bigger than l'.'O the pulsar light is some- 
what spread into the ring channels. Evidence of this effect 
is seen in the light curves when the intensities measured by 
individual ring fibres are folded with the pulsar phase and 
pulsed emission is detected. Additionally, there may be con- 
tributions of photons coming from the outer knot (not very 
clearly seen in the Fig.[l}, but an indication of it can be seen 
on the south-east side in between the fibres 3 and 4. It is 
also noteworthy to mention that the ring fibres see different 
patches of nebulosity. One third of all counts are recorded 
in the central fibre. The background fibres contribute to the 
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total number of counts in the percentage of 9, 10, 14, 11, 12, 
11 for channels from 1 to 6, respectively. 

To obtain the polarisation characteristics of the pulsar 
neighbourhood we assumed that within the DC phase range 
the contribution of the pulsar emission to the ring fibres 
is minimal. The integrated pulsed contribution in each fibre 
(spill-out from the pulsar) with respect to the total counts in 
the fibre is on the level of 3.1%, 4.5%, 4.4%, 2.5%, 2%, 1.9% 
for the channels from 1 to 6, respectively. Therefore, dur- 
ing our background (nebula) calculations we consider only 
light coming within the 'off-pulse' phase range, i.e. 7% of 
the whole rotational cycle of the Crab pulsar. Obtained po- 
larisation degree and position angles for each of the single 
OPTIMA apertures are given in Tab. [7] as well as illustrated 
in the Fig. [3] We compared our results with previous ones by 
over plotting them on the polarisat ion sky map of the very 
close neighbourhood presented by ISmith et al.l ll 19881 ). By 
averaging the Stokes parameters over all background chan- 
nels we get p = 9.71% ± 0.08% and 6 = 139.8° ± 0.2° for 
the region surrounding the pulsar. Close to the pulsar the 
nebular polarisation is quite uniform (~ 9 — 11%) but the 
position angles change steadily with radial distance. Two to 
three arc seconds from the pulsar the mean value is around 
140° but beyond five arc seconds the p osition angle exceed s 
155° and it is very position-dependent (|McLean et al.lfl983T l. 

5.2 Polarisation characteristics of the Crab pulsar 

The Crab pulsar is detected at all phases of rotation, i.e. 
also in the so-called 'off-pulse' phase with an intensity 
of about 2% compared to the maximum intensity of the 
main pulse. As mentioned before, the measured level of 
the DC differs from author to aut hor. For examp l e ver 
early measurements performed by IPeterson et al.l (|l97 
gives the 'unpulsed background' from the Crab pulsar on 
the level of 3.6% of the main peak i n tensit y. Much lower 
values were obtain ed by Ijones et alJ ()l98ll ) and later by 
ISmith et all (Il988fl : 0.6% and 1.2%, respectively. In addi- 
tion IPercival et all (|l993h claims that the 'off-pulsed' flux 
has an intensity less than 0.9% of the peak flux (visible and 
UV data from HST, 2a upper lim it), whereas a fractional 
flux derived bv lGolden et~aLl l|2000h from photometric anal- 
ysis gives ~ 1%. 

For the background (nebula and sky) subtraction we 
took the averaged Stokes parameters (Idc,Qdc,Udc) 
recorded in the ring fibres over the 'off-pulse' phase. The 
pulsar Stokes parameters (/, Q, U; Fig. IA1[) as a function of 
its rotational phase are derived after subtracting from the 
central channel measurements the steady nebular compo- 
nent. The colour coded Stokes parameters Q, U as a vector 
diagram are shown in the Fig. U Colours refer to the pulse 
phases as indicated in the input box and the scale is such 
that I — 100 at the maximum light, i.e maximum of the 
MP. The MP and IP maxima are indicated with black open 
squares. Points belonging to the MP phases follow an outer 
ellipse (upper panel), whereas these belonging to the IP an 
inner one (bottom panel). In both cases the direction of 
increasing pulsar phase is counter-clockwise. Noteworthy is 
that already from the Stokes parameters one can see that 
there is sudden change in the pattern near the radio phase, 
i.e. where the red points change to the black ones. As the 
next step of the data analysis the polarisation characteris- 



tics, the position angle and the degree of polarisation of the 
E- vector, are calculated from the Stokes parameters (Q. IA2|) . 
Results plotted with a different time resolution are presented 
in Fig. g] and Fig. [5] 

Optical emission from the Crab pulsar is highly po- 
larised, especially in the bridge and the 'off-pulse' phases. 
The position angle and the degree of linear polarisation as 
functions of rotation phase show well determined properties: 

a) The polarisation characteristics of both, the MP and 
the IP components, are quite similar. 

b) The polarisation degree reaches a minimum at phase 
close to the radio main peak; the minimum is not aligned 
with the optical peak. 

c) There is a well defined bump in the polarisation degree 
on the rising flank of the MP. 

d) There is an indication of such a bump also for the IP 
(especially after DC subtraction). 

e) The position angle swings through a large angle in both 
peaks: after subtraction of the DC component the angle 
swing is 130° and 100° for MP and IP, respectively. 

f) The position angle at the bridge and 'off-pulse' phases 
is constant. 

g) The position angle slope changes dramatically at phases 
0.993 (MP maximum) and 1.0 (radio peak). 

h) The trailing wing of the MP (phase range 1.0-1.03) 
shows a linearly increasing degree of polarisation. This fea- 
ture turns into a bump shape after DC subtraction. There 
is a slight indication of the same behaviour for the trailing 
wing of the IP. 



5.3 Polarisation characteristics of the Crab pulsar 
after DC subtraction 

The apparent constancy of the position angle within the 
phase range 0.78-0.84 (Fig. 2] right upper panel) may sug- 
gest that the optical emission from the Crab pulsar consists 
of two components - pulsed and unpulsed. The pulsed com- 
ponent is characterised by a highly variable position angle 
and polarisation degree. The unpulsed (i.e. DC) component 
is characterised by constant intensity on the level of 2% of 
the main pulse intensity, fixed 9 ~ 119°, and a degree of po- 
larisation on the level of 33%. The source of emission of the 
DC is unknown. There are various ideas from where in the 
magnetosphere or the nebula this component might arise. It 
is also possible that the inner knot (located only 0'.'65 apart 
from the pulsar and being a pers istent feature throu ghout 
the sequences of the HST images, iHester et al. Irt2002rn con- 
tributes to the 'off-pulse' emission, although Golden et al.l 
(2000) claim that the Crab image during the off-phase is 
compatible with an unresolved point source. Assuming that 
the unpulsed component is present at all phase angles and 
has constant polarisation we obtained the polarisation char- 
acteristics of the 'pulsed component' separately by subtract- 
ing the respective Stokes parameters Idc, Qdc, Udc- Re- 
sults are presented in the Fig. [6] and [7] After subtracting 
the 'unpulsed component' the position angle and polarisa- 
tion degree in the phases where the intensity is very low are 
not well defined. The values of 6 and p become very noisy 
because the Stokes parameters go basically to zero for these 
rotational phases. 

The polarisation degree of the 'off-pulse' component ob- 
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Figure 4. Left column: Stokes parameters Q, U as a vector diagram. Colours refer to the pulse phases as indicated in the input box 
of upper panel. The scale is such that I = 100 at the maximum light. Phases of the MP and IP maxima are indicated with black open 
squares. Points belonging to the MP and IP follow in a counter-clockwise direction on outer and inner ellipses, respectively. Bottom 
panel shows a zoomed region around the IP phases. The dashed lines indicate the origin of the scaled Q,U coordinates. Any error in 
background polarisation would correspond to a shift in the origin in this plot. The lower plot shows an zoomed region of the upper panel. 
Right column: The Crab pulsar position angle as a function of rotational phase (upper panel). Changes of 8 are aligned with the MP 
maximum of the optical light (vertical dot-dashed line), but also with the zero phase, i.e. with the radio phase of the MP (vertical solid 
line, see Fig. [5] for details). The Crab pulsar polarisation degree as a function of rotational phase (bottom panel). Minimum of p is for 
the radio phase (vertical solid line), and not for the maximum phase of the optical MP (vertical dot-dashed line, see Fig. [5] for details). 
For both panels a bit more than one rotation is shown for clarity. The pulse profile (solid line) and DC phase range (dashed region) are 
indicated also. The DC component in the phase range 0.78-0.84 has P. A. ~ 119°, and P.D. ~ 33%. One sigma errors are show. 



tained in thi s work differs from the value s presented by othe r 
authors: e.g. lJones etHI (|l98ll ): ~ 70%; ISmith et al.l i|l98Sl ): 
~ 47% ± 10%. This might be caused by two reasons. Firstly, 
different groups treat background subtraction in different 
ways. Additionally, all mentioned observations were taken 
during different epochs and with different instrumentation. 
Secondly, the observed variation might be caused by the 
intrinsic mechanism of the pulsar and/or nebula radiation. 
It should also be noted that the predicted d ecrease of th e 
optical luminosity is ~ 1% each two years l|Pacini|[l97ll ), 
therefore one wo uld expect a r educt ion of 0.13 mag over the 
22 years between I Jones et al.l (|l98ll ) and our observations. 



5.4 Alignment between optical and radio 
wavelengths 

Precise timing of pulsar light curves throughout the elec- 
tromagnetic spectrum can be used to constrain theories of 
the spatial distribution of various emission regions and their 
specific propagation delays. In the radio regime this concept 
has been applied in the so called frequency mapping analy- 
sis. From recent observations at different energies, it became 
clear that the Crab pulsar emission maxima of MP and IP 
are not aligned in phase at different wavelengths, from ra- 
dio to the 7 energy range. Com paring the visible an d UV 
light curves obtained from HST, iPercival et all (|l993l ) were 
among the first to show that the phase separation between 
the two peaks, ipso facto the phases of the pulse maxima, 
change with energy. Since then many authors have mea- 
sured this effect in an attempt to understand its relation 
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Figure 5. Polarization characteristics of the Crab pulsar - P.A. (top row) and P.D. (bottom row) - as a function of rotational phase. 
Data binning is 1000 and 3000 bins per cycle (corresponding to 33.5 fis and 11 /ts per bin) for odd and even columns, respectively. The 
first two columns are zoomed around the MP phase (using different scales), whereas the next two are zoomed around the IP phase. 
Dot-dashed vertical line indicates the optical MP maximum phase (=0.993), as well as the IP maximum phase (= 1.398). Solid vertical 
line indicates the radio peak phase (=1.0). For clarity the optical light curve of the Crab pulsar is over plotted (solid line). 
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Figure 6. The same as in Fig.[5]but after DC subtraction. 



to and impact on the emission mechanism. However, the 
techniques used by different authors to measure the phase 
separation (i.e. the phases of the peak maxima) have var- 
ied. Therefore, this might cause method-dependent biases. 
lEikenberrv fc Fazid (|l997h showed that the peak-to-peak 
separation appears to be a more or less smooth function of 



energy from infrared to 7-ray energies. The separation de- 
creases from 0.4087±0.0003 to 0.398±0.003 with energy over 
the range from 0.5 eV to 1 MeV, respectively. There is some 
evidence of a turnover or a break in this trend at energies of 
0.7 eV (H band pass filter). No default method exists for de- 
termining the position of the peak of the profile. For our pur- 
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Figure 7. The same as Fig. [4] but after DC subtraction. 



pose and using our high statistics light curve it was enough 
to determine the peak phases just by looking for the max- 
imum intensity. This gives us the phase values for the MP: 
0.993 ± 0.001 and the IP: 0.398 ± 0.001 . Our peak-to-peak 
separation is on the level of 0.4050=1=0.00 14. It is in very good 
agreem ent with the values obtained by lEikenberrv fc Faziol 
1 19971 ) for the visual pass band, i.e. 0.4057 ± 0.0003, and 
from a previous OP TIMA observation of 0.4060 ± 0.0003 
jStraubmeierl (120011 1). 

Simultaneously with the optical observations we per- 
formed radio observations at Jodrell Bank observatory, 
which allows us to compare these wavelength ranges. Pulsed 
radio emission at 610 and 1400 MHz is detectable only 
around the peaks of emission due to the strong background 
from the nebula. Fig.[8]shows the radio light curves. At 610 
MHz the MP is preceded by the so-called precursor peak, 
which is mostly visible at lower frequencies and is not de- 
tected at 1400 MHz. The nature of this radio precursor is 
unclear - some researchers proposed that it is the proper 
emission from the polar gap of the pulsar while the main 
peaks are ge nerated higher up in the magnetosphere (e.g. 
lRankin|[l990l) . Recently, Petrova (|2008h showed that differ- 
ent components in the Crab pulse profile may be induced 
by the scattering from different harmonics of the particle 
gyrofrequency that takes place at different magnetospheric 



altitudes. This and the rotational effect give rise to the 
components. In this model, the low frequency component 
is formed by the scattering from the first harmonic of the 
gyrofrequency into the state below the resonance, whereas 
the precursor by the scattering between the states below 
the resonance. This model is well supported by the radio 
polarization data. 

As can be seen in Fig. [8] the optical maxima of the MP 
are leading the radio pulse. A similar lead has also been 
found at X- and 7-ray energies. T he reported value s of the 
time lag are as follo w: 344 ± 40 us dRots et alj|200"j . RXTE 
data), 280 ± 40 us dKuiper et alll2003l . INTEGRAL data), 
and 241 ± 29 us (|Kuiper et alj 120031 EGRET data). The 
uncertainty of the latter value does not include the EGRET 
absolute timing uncertainty of better than 100 us. At opti- 
cal wavelengths the situation is quite different and does not 
provide such a coherent picture. The quantitative amount 
of the lag between optical and radio peaks, and whether 
it exists at all and at all times, is controversial when we 
compare results obtained in different investigations. Several 
authors reported the optical p eak leading t he radio peak 
by a time shift of 1 40 ± 78 us (|Sanwallll999l ). 100 ± 20 us 
dShearer et al.l l2003h. an d 273 ± 100 us rece ntly found by 
lOosterbroek et al. ( 2006). On the other hand iGolden et al.l 
(2000) reported that the optical pulse trails the radio pulse 
by about 80 ±60 us. Additionally. iRomani et all (120011 ) con- 
cluded that both, radio and optical, peaks are coincident to 
better than 30 us, but this analysis did not take into account 
the uncertainty of the radio ephemeris being on the level of 
150 us in the error calculations. 

From our measurements we conclude that the optical 
main peak is leading the radio peak by a time shift of 
231 ± 68 us. The uncertainty in this value is composed from 
an uncertainty of 33 us in the determination of the opti- 
cal peak of the MP and from 60 us in the radio ephemeris. 
Our value of the optical phase difference between the MP 
and IP of 0.4050 ± 0.0014 is cons istent with the late s t opt i- 
cal measurements carried out by lOosterbroek et alj (2006), 
who obtained 0. 4054±0.0004. Both v alues , as it has already 
been shown by Ei kenberrv fc Faziol (|l997n . are not consis- 
te nt with the X-ray results, e.g. obtained from RXTE data 
bv lRots et all j|2004h of 0.4001±0.0002. This implies that the 
details of the pulse profile in X-rays and in the optical do- 
main are different. In a simple geometrical model (ignoring 
relativistic effects) a time shift of ~ 231 us indicates that 
possibly the optical radiation is formed ~ 70 km higher 
in the magnetosphere than the radio emission. The differ- 
ence in phase of 0.007 could also be interpreted as an angle 
between the radio and optical beam of ~ 2.5° (neglecting 
aberration and magnetic sweep back). 

Another striking correlation between the radio intensity 
profile and the optical polarisation can be seen in Fig. [8] the 
radio precursor seems to be perfectly aligned with the bump 
in the degree of optical polarisation. During this phase of the 
leading wing of the optical MP the position angle change is 
also characterised by a nearly linear swing. At the present 
stage of modelling the coherent and incoherent emissions 
from a pulsar magnetosphere, where both processes are gen- 
erally treated independent of each other and mutual inter- 
actions have not been investigated in depth, it is premature 
to speculate on the origin of these observational results. 
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Figure 8. Optical polarisation characteristics of the Crab pulsar (P. A. - top row, P.D. - bottom row) compared with the pulsar radio 
profiles. Radio profiles obtained at radio frequencies of 1400 MHz and 610 MHz are shown as a solid line and dashed line, respectively. 
Left column shows zoom around the main pulse phase, whereas right column around the inter pulse phase. Points indicate the optical 
polarisation measurements, while the dotted line shows the optical intensity profile. 



Table 8. Polarisation degree (P.D.) and position angle (P.A.) at 
three selected Crab pulsar phases, i.e. the phases of both optical 
maxima, IP (0.398) and MP (0.993), as well as, the radio phase 
(0.000). Values before and after DC component subtraction are 
shown. 









After DC Subtraction 


Phase 


P.D. [%} 


P.A. [°] 


P.D. [%] 


P.A. [°] 


0.398 


7.6 ±0.3 


109.9 ± 1.1 


5.8 ±0.3 


106.0 ± 1.6 


0.993 


5.9 ±0.2 


119.8 ±0.8 


5.3 ±0.2 


119.9 ±0.8 


0.000 


3.3 ±0.2 


149.6 ± 1.6 


3.0 ±0.2 


157.5 ± 1.7 



6 SUMMARY AND DISCUSSION 

The Crab pulsar emits highly anisotropic radiation which 
spans a wide range of wavelengths, from radio to extreme 
7-rays. Knowledge of polarisation characteristics of this ra- 
diation is of fundamental importance in our attempts to 
find the mechanisms responsible for Crab's magnetospheric 
activity. Good quality X-ray and 7-ray polarimetry with 



satellite observatories is expected to be available for pulsar 
studies (among other types of objects) in the near future. 
Present-day state of instrumentation allows to carry out op- 
tical polarimetry of this object with unprecedented quality. 
Our project to study the Crab pulsar with OPTIMA at NOT 
is, to the best of our knowledge, the most recent and most 
complete one. The observations with 11 [is time resolution 
are an order of magnitude better than the previous best 
observations. We have completely resolved the polarisation 
characteristics of both peaks of the Crab pulsar, MP and IP, 
in the optical pass bands (see Fig.[S}. Moreover, we were able 
to better characterise the polarised emission between the 
peaks, i.e. the bridge as well as the DC ('off-pulse') region. 
We find that the MP of the Crab pulsar arrives 231 ± 68 /is 
before the peak of the radio pulse. 

The phase averaged polarisation degree of the Crab 
pulsar amounts to 9.8% with a position angle of 109.5°. 
After the DC subtraction it is 5.4% and 96.4°, respec- 
tively. Minimum polarisation degree occurs at the phase 
of 0.999, very close to the radio pulsar phase, where p = 
3.3% ± 0.2%, = 149.6° ± 1.6° before the DC subtraction, 
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and p = 3.0% ± 0.2%, 6 = 157.5° ± 1.7° after DC subtrac- 
tion. During the IP the minimum value of p is on the level of 
5%, before DC subtraction. Similar to the MP case, it is also 
shifted (with respect to the phase of the optical maximum) 
to the leading wing of the pulse. The value of p changes af- 
ter DC subtraction to ~ 2%. This means that the situation 
is inverted. Before the DC subtraction the minimum of po- 
larisation degree is reached during the main pulse, whereas 
after this subtraction it is observed during the inter pulse. 

The Crab is the pre-eminent example of a high-energy 
source with polarised optical emission - if the pulsar is not 
phase resolved, the polarisations still averages out to about 
10%. Other optical pulsars have been measured to have an 
average polarisation also up to 10%. Polarization can be 
generally expected if the radiating particles are confined, e.g. 
by magnetic forces, to anisotropic distributions. Detection 
of polarised objects in the field of unidentified 7-ray sources 
could be therefore a valuable tracer to aid identification. 

Our result s agree generally well with previous mea- 
surements fe.g. lSmith" et al.ll 19881 ; iKanbach etal]|2005l ) . but 
they show details with much better definition and statistics. 
The behaviour of 9 as a function of phase observed for the 
Crab pulsar at optical wavelengths (Fig. [4| differ from those 
observed at radio wavelengths (e.g. Mo ffett fc Hankinsl 19991 ; 



iKarastergiou et al.ll200i ; ISlowikowska et al.ll2005h . Two fac- 
tors may be responsible for this difference: different propa- 
gation effects, and different intrinsic emission mechanisms. 
In particular, the former factor plays an essential role in 
such high-energy emission models like the outer gap model 
or the two-pole caustic model: in both models high-energy 
emission comes from a very wide range of altitudes, contrary 
to radio emission which originates within a narrow range of 
altitudes. For comparison the light curves and polarisation 
characteristics obtained within the framework of three high 
energy magnetospheric emission models of pulsars, i.e. the 
polar cap model, the two-pole ca ustic model, and t he outer 
gap model are shown in Fig. [TO] (|Pvks et al.ll2004alE). 

The two-pole caustic model l|Dvks fc Rudakll2003l ) pre- 
dicts fast swings of the position angle and minima in the 
polarisation degree, similar to what is observed. Polar cap 
model and outer gap model are not able to reproduce the 
observational polarisation characteristics of the Crab pulsar. 
Another model, placing the origin of the pulsed optical emis- 
sion from th e Crab in a striped pulsar wind zone has been 
proposed by IPetri fc Kirkl (I2005T ). This model features also 
polarisation characteristics that bear a certain resemblance 
to the observ ations. 

Recently. iTakata et al. (2007) attempted to simultane- 
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Figure 10. The optical light curve, the position angle and the degree of polarisation calculated with the following m odels of high energ y 
radiation from pulsars, from left to right: the polar cap mod el, the two-pole caust ic model, and the outer gap model l|Dvks et al.ll20041 
In the first column, the data from Calar Alto observations (Kanbach ct al. 2005) are shown for comparison. Courtesy J. Dyks. 



ously reproduce all known high energy emission properties 
of the Crab, including its optical polarisation characteris- 
tics within the framework of a modified outer gap model. 
The model is restricted to synchrotron emission due to sec- 
ondary and tertiary electron-positron pairs which are ex- 
pected in different spatial locations of the 3D gap; inter- 
nal polarisation characteristics are calculated with particu- 
lar care. Yet, the calculated polarisation properties for op- 
tical light hardly reproduce the observed properties. How- 
ever, similarly like in the case of the two-pole caustic model 
(Dyks et al. 2004a bL the p olarisation characteristics ob- 
tained in Taka ta et al. | i|2007f ) become more consistent with 
the Crab optical data after the DC component is subtracted 
(as men tioned in the previo us section). An important out- 
come of I Takata et alj (|2007l ) is that it offers the energy de- 
pendence of the polarisation features, covering the energy 
range between 1 eV and 10 keV. 

The models of pulsar magnetospheric activity are based 
on various (sometimes ad hoc) assumptions and different 
boundary conditions. Those lead to the model differences 
in the macro scale (spatial extent of accelerators and emit- 
ting regions), as well as in the micro scale (specific radia- 
tive processes). The former include polar gaps, slot gaps, 
caustic gaps, outer gaps and striped winds. The latter in- 
clude e.g. curvature radiation, synchrotron radiation and 
inverse Compton scattering. In consequence, the models dif- 
fer significantly in the resulting 'observed' radiation prop- 
erties: light curves, energy spectra, and - last but not 



least - polarisation dDvks et alj|2004allbt iPetri fc Kirkll2005l : 
iRomani fc Yadigaroglulll995l ; iTakata et al.ll2007l ). 

Linear polarisation characteristics in the high energy 
domain (optical, X-rays and 7-rays) is considered as a pow- 
erful tool which may lead to a breakthrough in our un- 
derstanding of pulsar emission mechanism. For this reason, 
the optical high time resolved polarisation properties ob- 
tained for the Crab pulsar have attracted particular atten- 
tion due to their un iqueness. Some mo dels, like the two- 
pole c austic model (|Dvks et all 2004a|3), the outer gap 
model (jRomani fc Yadigaroghl 19951: ITakata et al.|[2007l ) or 
the striped pulsar wind model ( Petri fc Kirkll2005h are able 
to reproduce (very roughly) some of these properties, e.g. 
(a) and (e). However, a fully convincing explanation of the 
properties listed from (a) to (h) is beyond the reach of all 
above-mentioned models. 
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APPENDIX A: 

Here, we present step-by-step our polarisation data analysis based on the method described bv lSparks fc Axo n (1999). The 
rotating polarisation filter (RPF) inside the OPTIMA instrument provides polarimetric data that represent a series of 'images' 
of an object taken through 180 sets of linear polarisers, when we bin the continuous rotation of the RPF into discrete one 
degree intervals. A single polariser is not a 100% perfect polariser, but its characteristics are well established, and this is 
essential for the chosen data analysis method. From an input data set of 180 idependent intensities / (measured in counts) 
and their errors (07 = \/7), corresponding to a set of observa tions through 18 ident ical but not perfect polarisers we derive 
the Stokes parameters, following the case of n polarisers after ISparks fc Axon! l|l999h . 

Linearly polarised light requires measurement of three quantities to be fully characterised. There are various ways of 
expressing this. The most common one involves the total intensity of the light I, the degree of polarisation p, and the position 
angle 9. An intermediate stage between the input data and the solution of polarisation quantities are the Stokes parameters 
(I, Q, U) that are related through: 

Q = Ip cos 26i, U = Ip sin 29, ( Al) 

or equivalently, 

(Q2 + U2) 1/2 a 1 (U\ . 

p =y^ — _) — _ = _ arctan ( \ (A2 ) 

These quantities describe all intrinsic properties of linearly polarised radiation. We would like to underline that they should 
not be confused with the properties of the polarising elements of the polarimeter, which in our case are the properties of 
the measured intensities in each one degree intervals of the RPF, i.e. in each of 180 polarisers. These three quantities that 
characterise fully the behaviour or response of the linearly polarising element are: 

- its overall throughput (hereafter t), in particular to unpolarised light, 

- its efficiency as a polariser (hereafter e), i.e. the ability to reject and accept polarised light of perpendicular and parallel 
orientations, 

- the position angle of the polariser (hereafter (j>). 

There are a vari ety of conventions co mmonly used to present these quantities dMazzuca et al.lll998t ). Here, we adopt the 
convention after ISparks fc Axon|[l999l . i.e. the output intensity of a beam with input Stokes parameters (I, Q, U) passing 
through a polarising element is given by 

h = i*fc[/ + e fc (cos2<^Q + sin2<^£/)] (A3) 

where tk is related to the throughput to unpolarised light, e k is the efficiency of the polariser, and 4>k is the position angle 
of the polariser k. During our measurements we always use the same Polaroid, just the position angle of the filter changes, 
therefore = e = 0.998, t& = t = 0.32 (reference: Polarisation Filter Type VIS 4 K, Linos Photonics), and <f)k takes values 

from to 179 with 1 degree steps. 

Following the equations given bv ISparks fc Axor] Jl99sh we define the new three-component vector for the effective 
measurements: 

T ll V - " • Iktk r ll \^ IkXk^k COS 2(j> k T ll \^ Iktk^k SU1 2<f>k .x 

h -^^' h — — ' — OV 2 — • (A4) 

The error estimate for each measured intensity (^counts) is based on Gaussian statistics; therefore Ik = Ok" and the 
corresponding factors reduce to unity. The results of calculating the effective intensity components I\,I'<1 and I3 are shown 
in Fig. IA2I In our case these values are constant. 7" is just a sum over the same 180 polarisers with the same throughput 
to unpolarised light, i.e. t = 0.32. Therefore we get: l'{ = 180 x 0.32 = 57.6. Whereas, I' 2 ' = 180te £ cos 2c/> fc and = 
180te £ sin 20fe , where £ denotes a sum over index k of the 180 polarisers. Knowing that e = 0.998 and the integrals of 
cos20fe and sm2(j>k over the range from 0° to 180° we expect 1% ~ /" X] cos 2</>fc and J3 ~ I" X shr 2</>fc, within available 
numerical precision, to be very close to zero. 

Similarly, we can define a three-component vector of effective transmittances 

,// sr^ t k 2 ,;, «— t k 2f ^k cos2(f>k .a v-^ t k 2 £k sin 2</> fe 

ti = y. — j ' *2 = > - 2 , h = > , (A5) 

and the vector of effective efficiencies 



E*fc 2 / CT fe 2 V V o"fc 2 / V o-fc 

1 



E*fc 2<E fc cos2 0fe/o- fc 2 V V o-fc 2 / V °k 



^2 ~^2 e k 2 cos 2 2<f> k J + ~^ e k 2 sin 2</> fc cos 2<f> k 



2 



//^*fcl £fc2 8^2^ 008 2^) + (j2^k 2 sm 2 2cj,k 
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Figure Al. Stokes parameters /, Q, U derived from 180 polarisers h as a function of the Crab pulsar rotation phase (1 cycle = 500 bins). 
Normalisation was such that / = 100 at the maximum, and Q and U are normalised to J. These values corresponds to the values shown 
on the QU plane in the Fig. [4] 



as well as a vector of effective position angles 



1 arctan ( ^re fe sin 20 fc / ^rtfe cos 20 

2 o k l / < a k 2 

~ arctan ( ^re fc 2 sin 2(j) k cos 20 fc / Y^ ^^e k 2 cos 2 2<f> k 
2 \ ^— ' a k 2 / ' a k 2 

1 arctan ( ^-^e k 2 sin 2 2<f> k / Y^ ^-e t 2 sin 20 fe cos 20 fe 

2 <7fe 2 / ^ (Jfe 2 



The effective intensity, transmittances, efficiencies, and position angles are shown in Fig. IA2I By making these substitu- 
tions, the solution for the Stokes vector is given by 

(I,Q,U) = B[ tf/(0.5tf) (A6) 



V tf/(0.5#) / 

/ e 2 'e 3 'sm(20 3 '-20' 2 ') ei'e^'sin(2^i'-2^') 44 sin(20 2 ' - 20'/) \ 
B=i £2 sin 202 - £3 sin 203 e 3 sin 20 3 - e" sin 20'/ 4 sin 20" - 4 sin 20 2 ' (A7) 
\ e 3 ' cos 20 3 ' - 4 cos 20 2 ' e'i cos 20'/ - 4 cos 20 3 ' e 2 ' cos 20 2 - 4 cos 20i' / 

= 44 sin(202 - 20") + 44 sin(203 - 20 2 ) + 44 sin(20" - 203). (A8) 

can be used immediately. In this way we obtained the Stokes parameters shown in Fig. IA1I 

W e calculated the covariance matrix (Fig. I A3f) defined as the inverse of matrix C given by Eq. 8 in Sp arks fc Axon! 
(1999). Following the error propagation equation we calculate the uncertainties of the polarisation degree and the position 
angle according to: 



where 



and 



< - + W'°t, + + 2 ^)(^) + + (A9) 

^-^(f) 2 + 4(g) 2 + 2^(§)(g). (A10) 
where appropriate standard deviations are the components of the covariance matrix. In general the diagonal terms 
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dominate the uncertainties. Covariant terms make the maximum contribution on the level of 5% and 1% to the a p and 
respectively. 
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